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We study ultracold fermionic atoms trapped in an optical lattice with harmonic confinement 
by dynamical mean-field approximation. It is demonstrated that a supersolid state, where an 
s-wave superfluid coexists with a density- wave state with a checkerboard pattern, is stabilized 
by attractive onsite interactions on a square lattice. Our new finding here is that a confining 
potential plays an invaluable role in stabilizing the supersolid state. We establish a rich phase 
diagram at low temperatures, which clearly shows how an insulator, a density wave and a 
superfluid compete with each other to produce an interesting domain structure. Our results 
shed light on the possibility of the supersolid state in fermionic optical lattice systems. 
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Ultracold atomic gases have attracted much interest^ 
since the successful realization of Bose-Einstein conden- 
sation in a bosonic ^^Rb system.^ Optical lattices, formed 
by loading ultracold atoms in a periodic potential, have 
been providing an ideal stage for experimental and the- 
oretical studies of fundamental problems in condensed 
matter physics. ^^ Owing to its high controllability in 
interaction strength, particle number, and other param- 
eters, many remarkable phenomena have been observed 
such as the phase transition between a Mott insulator 
and a superfluid in bosonic systems.'' More recently, a 
fermionic gas in the optical lattice has been a topic of 
extensive study, which has successfully lead to the obser- 
vation of superfluidity in the case of attractive interac- 
tions.^ 

One of the interesting questions for such a fermionic 
optical lattice is how an 5-wave superfluid (SSF) state 
coexists or competes with a density wave (DW) state, 
where the latter can be regarded as a sort of solid state. 
This provides an important issue in condensed matter 
physics, since it is directly related to a hot topic of cur- 
rent interest, the so-called supersolid state. The existence 
of the supersolid state was suggested for ^He experimen- 
tally,^ and this pioneering work has stimulated theoret- 
ical investigations on bosonic systems ^^^^ and Bose- 
Fermi mixtures. ^^ As for fermionic systems, the optical 
lattice can be a potential candidate for it. However, it has 
not been clarified how the supersolid is stabilized in the 
optical lattice except for the one-dimensional system^^' ^^ 
although the existence of the SSF and DW states has 
been discussed. ^^^^ 

According to previous studies of the attractive Hub- 
bard model on a periodic lattice without a confining po- 
tential, ^^^^ in systems on bipartite lattices, except for 
a one-dimensional case, the DW and SSF ground states 
are degenerate at half filling, which means that the su- 
persolid state might be realizable in principle. However, 
the degenerate ground state is unstable against pertur- 
bations. For example, away from half filling, the super- 
solid state immediately changes to a genuine SSF state, 
where a BCS-BEC crossover has been discussed.^^'^^'^^ 



In the optical lattices, we have an additional confining 
potential, which makes the situation different from the 
homogeneous bulk system. This naturally motivates us 
to address the question whether the supersolid state can 
be realized in the optical lattice. 

In this study, we demonstrate that the supersolid state 
can indeed be realized in a fermionic optical lattice with 
attractive interactions. In particular, it is found that a 
confining potential plays an important role in stabilizing 
the supersolid state; it makes the supersolid state robust 
against perturbations in contrast to that in homogeneous 
systems. This suggests that the fermionic optical lattice 
can be a potential candidate for the realization of the 
supersolid state. 

We consider ultracold fermionic atoms, which may be 
described by the Hubbard model with confinement as 



H 






c]^c^ 



uE^ 



ni^Tiii + 



Vb^iifni,, (1) 



where Cicr{c\^) annihilates (creates) a fermion at the zth 
site with spin a and n^cr = ScrQa- ^ij is the nearest- 
neighbor hopping, U the attractive interaction, and Vb 
the curvature of a harmonic potential. Here, Ri is the 
distance measured from the center of the system. 

The ground-state properties of the Hubbard model on 
inhomogeneous lattices have theoretically been studied 
by various methods such as Bogoljubov-de Gennes equa- 
tions,^^ Gutzwiller approximation,^^' ^^ and variational 
Monte Carlo simulations.^^ Although ordered states are 
described properly in these approaches, it may be diffi- 
cult to deal with the coexisting phase like a supersolid in 
the strong correlation regime. The density matrix renor- 
malization group method^^' ^^ and quantum Monte Carlo 
method^^'^^ are efficient for one-dimensional systems, 
but may be difficult to apply to higher-dimensional sys- 
tems. We here use dynamical mean- field approximation 
(DMFA) ^^^^^ which incorporates local particle corre- 
lations precisely, thus enabling us to obtain reliable re- 
sults if spatially extended correlations are negligible. In 
fact, the method has successfully been applied to some 
inhomogeneous correlated systems such as the surface^^ 
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or interface of Mott insulators^^ and repulsive fermionic 
atoms. ^^'^^ An advantage of this method is to treat the 
SSF and DW states on an equal footing in the strong 
correlation regime. 

In the framework of DMFA, the lattice Green's func- 
tion is described in terms of the site-diagonal self-energy 
ti {iiJn) as 

where cFz is the z component of the Pauli matrix, cfq 
the identity matrix, /i the chemical potential, uon = 
(2n + l)7rT the Matsubara frequency, and T the tem- 
perature. A DMFA self-consistent loop of calculations is 
iterated under the condition that the site-diagonal com- 
ponent of the lattice Green's function is equal to the 
local Green's function obtained from the effective im- 
purity model. ^^^^ When DMFA is applied to our inho- 
mogeneous system, it is necessary to solve the effective 
impurity models L times by iteration, where L is the sys- 
tem size. For this purpose, we use a two-site approxima- 
tion.^^' ^^ Although the effective bath is replaced by only 
one site in the method, it has the advantage in taking 
into account both low- and high-energy properties rea- 
sonably well within restricted numerical resources. ^^'^^ 

In the following, we consider the square lattice with 
harmonic confinement as a simple model for the super- 
solid. We set t as a unit of energy, and fix the curva- 
ture of the potential and the total number of atoms 
as Vo = 0.023 and N - 300(7V^ - 150). In our sys- 
tem, the distribution of particles and order parameters 
is spatially modulated, which is optimized in the frame- 
work of DMFA. We thus calculate the density profile 
(^icr) = 2T^^^QRe[Gicr {i^n)] + ^ and the distribution 
of the pair potential A^ = 2T^^^QRe[F^ (iujn)], where 
Gia{iuJn)[Fi{iuJn)] is the normal (anomalous) Green's 
function for the ith site. Note that A^ represents the 
order parameter for the SSF state. 

The obtained results at T = 0.05 are shown in Fig. 
1. In the weak coupling case {U = 1), fermionic atoms 
are smoothly distributed up to i? ~ 13, where R is the 
distance from the center of the harmonic potential. In 
this case, the pair potential is not yet developed [Figs. 1 
(a) and (b)]. Therefore, a normal fluid state with short- 
range pair correlations emerges in the region {R < 13). 
Increasing the attractive interaction, fermions tend to 
gather around the bottom of the harmonic potential, as 
seen from (n^cr) in Fig. 1 (c). Note here that the attrac- 
tive interaction causes an SSF state with finite A^ in the 
region with (n^cr) 7^ 0, as shown in Fig. 1 (d). This is con- 
sistent with the results obtained from the Bogoljubov-de 
Gennes equation. ^^ In this figure, we encounter a remark- 
able feature around the center of the harmonic potential 
(i? < 5): a checkerboard structure appears in the density 
profile (n^cr), while the SSF state is not suppressed com- 
pletely even in the presence of the DW state. This implies 
that the DW state coexists with the SSF state, i.e., a su- 
persolid state is stabilized in our optical lattice system. 
The characteristic properties of the supersolid state are 
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Fig. 1. (Color online) Density profile (riia) (left panels) and pair 
potential A^ (right panels) on square lattice at T = 0.05 when 
U = 1.0, 3.0, 5.0, 7.0 and 9.0 (from top to bottom). 



clearly seen in the case of /7 = 5, where the checkerboard 
structure appears in the doughnut-like region [Fig. 1 (e)]. 
Further increase in the interaction excludes the DW state 
out of the center. It is seen in Figs. 1(g) and l(i) that 
fermionic atoms are concentrated around the bottom of 
the potential for large U. In the region, two particles 
with opposite spins are strongly coupled by attractive 
interaction to form a hard-core boson, giving rise to an 
insulating state with (n^cr) ~ 1. We observe such behav- 
ior more clearly in Figs. 1(h) and l(j). 

To observe how the supersolid state is realized, we also 
show the spatial variations in (n^cr) and A^ in Fig. 2 as 
functions of R. It is found that (n^cr) and A^ describe 
smooth curves for i? < 3 and 10 < i? < 11, where the 
SSF state without the DW is realized. On the other hand, 
for 3 < i? < 10, two distinct amplitudes appear in (n^cr), 
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Fig. 2. (Color online) Open and solid circles represent {riia-) and 
Ai at U = 5, and solid squares represent (riia) in the noninter- 
acting case at T = 0.05. 



reflecting the fact that the DW state with two sublattices 
is reahzed. An important point is that the pair potential 
A^ is finite in the region although its profile is somewhat 
affected by the spatial variation in DW. We thus confirm 
that the supersolid state is realized in the doughnut-like 
region (3 < i? < 10). 

By performing similar calculations, we end up with the 
phase diagram shown in Fig. 3. There are several remark- 
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Fig. 3. (Color online) Phase diagram of attractive Hubbard 
model on optical lattice with Vq = 0.023 and A^ ~ 300. The 
density plot represents the amplitude of the s-wave pair poten- 
tial A^. The DW state is realized in the shaded area. The broken 
lines serve as a visual guide which distinguishes the region with 
a fractional particle density from the empty and fully occupied 
regions. 



able features in the phase diagram. First, we notice that 
the supersolid state, characterized by the coexistence of 
SSF and DW orders, is indeed stabilized in a finite re- 
gion ([/ :^ 3 ~ 7), which is surrounded by the pure SSF 
state extended in a wider region. In the pure SSF region 
without the DW order, we still observe an interesting 
behavior, i.e., a BCS-BEC crossover. When U is rather 
small, the weak attractive interaction stabilizes a BCS- 
type SSF state, where A^ is induced in the whole region 
with (riicr) 7^ 0. In this parameter region, the pair poten- 
tial is enhanced with an increase of the interaction U. 
On the other hand, in the strong coupling region, parti- 
cles form short-range pairing states. In fact, most of the 
particles condense around the center yielding the insulat- 



ing state, and the others form a EEC-type SSF state in 
the vicinity of R = 7. Further increase in the attractive 
interaction narrows the SSF region, and suppresses the 
amplitude of the pair potential. Therefore, the pair po- 
tential Ai has a maximum at approximately U = 12 and 
i? = 7, which may give a rough guide for the crossover 
region between the BCS-type and BEC-type SSF states. 
We finally note that our supersolid state found in the 
weak coupling region at approximately U = 5 is thus at- 
tributed to the coexisting state of the DW and BCS-type 
SSF states. 

Next, we would like to discuss finite-temperature prop- 
erties in more detail. Note that we employ mean- field 
approximation for ordered phases, so that the corre- 
sponding transition temperature is finite even in two 
dimensions. Nevertheless, some essential properties of 
the supersolid at finite temperatures can be captured 
by the preset treatment; the results may be applied to 
the case where a weak three dimensionality is introduced 
as should be in real experiments. Here, we focus on the 
case of /7 = 5 to clarify how robust the supersolid state is 
against thermal fluctuations. The DW state is character- 
ized by the checkerboard structure in the density profile 
(n^cr), so that Fourier transform nq[= ^{ui^) exp(zgi?i)] 
at g = (tt, tt) is an appropriate quantity to discuss its 
stability. In Fig. 4, we show n(7^^7^)/n(o,o) and the max- 
imum of the pair potential A^ax- Decreasing temper- 
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Fig. 4. (Color online) Maximum A^ and normalized ^(7r,7r)/''T'(o,o) 
when [/ = 5 as functions of temperature T. The inset shows 
{riicr) and A^ as functions oi R. Crosses, circles, triangles, squares 
represent the results at T = 1.0, 0.5, 0.2 and 0.05. 



atures, n(7^^7^)/n(o,o) becomes finite below Tdw{^ 0.5), 
where the DW state is realized. Once the DW is ordered, 
the corresponding spatial region with the checkerboard 
structure (3 < r < 10) is hardly affected by tempera- 
ture, as shown in the inset of Fig. 4. On the other hand, 
the maximum pair potential Amax starts to develop at 
Tssf{^ 0.2), which is lower than Tdw- Therefore, the 
supersolid state with finite rii^^^^^ and A^ax is stabilized 
below TssF' These results imply that the DW state is 
more stable than the SSF state against thermal fluctu- 
ations for these parameters. On the other hand, around 
the border in the supersolid region in Fig. 3, the DW 
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state becomes unstable. Therefore, to observe the super- 
sohd state experimentally, it may be necessary to find ap- 
propriate parameters that should stabilize the DW state, 
since the SSF state is rather stable in the wide parameter 
region. 

Finally, we emphasize the importance of the confining 
potential to realize the supersolid state in optical lattice 
systems. Let us recall that the DW state with a checker- 
board pattern emerges when a commensurability condi- 
tion, e.^., half-filling, is satisfied for the particle density 
at least locally. Therefore, in order to stabilize the su- 
persolid state, it is essential to form a domain where the 
commensurability condition is met approximately. In our 
results presented here, such a domain is indeed formed in 
the doughnut-like region (Fig. 1). If the strength of the 
confining harmonic potential is decreased with N/V ^ 1 
fixed, where V is an effective system size, the transition 
temperature Tdw approaches zero, and finally the super- 
solid state changes to a pure SSF state, since the domain 
that satisfies the half-filling condition disappears. There- 
fore, we claim that a confining potential, which gives rise 
to an inhomogeneous distribution of the particle density, 
plays a key role in stabilizing a supersolid state in the 
fermionic optical lattice. This in turn demonstrates that 
the optical lattice system could be a potential candidate 
for realizing a supersolid state. 

We wish to comment on accessible experimental 
parameters for observing the supersolid state. The 
depth and curvature of the lattice and harmonic po- 
tentials can be controlled by adjusting the inten- 
sity and frequency of lasers. The hopping integral 
t and the attractive interaction are then given as 
t/En - 4^-i/2(^o/^^)3/4g-2(^o/£;K)^/^ and U/Er - 
— {S/7ry^'^askL{vo/Eii)^^^^ where '^o and /cl are the in- 
tensity and wave number of the laser for the lattice, Er 
the recoil energy, and a^ the s-wave scattering length. ^^ 
In the paper, we have seen that the supersolid state ap- 
pears in the vicinity of the BCS-BEC crossover region, 
which implies that the supersolid state is experimentally 
in accessible regions. Therefore, the supersolid state is 
expected to be observed by tuning these experimental 
parameters in the near future. 

In summary, we have investigated the fermionic at- 
tractive Hubbard model in an optical lattice with har- 
monic confinement. Using DMFA, we have obtained a 
rich phase diagram on a square lattice, which has an 
interesting domain structure including the SSF state in 
the wide parameter region. In particular, we have found 
that the supersolid state, in which the SSF state coexists 
with the DW state, is stabilized at low temperatures. It 
has also been elucidated that a confining potential plays 
a key role in stabilizing the supersolid state. There are 
many interesting problems to be explored in this context. 
An imbalanced fermionic system with N^ ^ N^ may be 
particularly interesting, since the so-called Fulde-Ferrell- 
Larkin-Ovchinnikov superfiuid state with a periodically 
modulated order parameter might emerge and compete 
with the DW state, giving rise to a novel supersolid state. 
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